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Abstract 

We calculate DIS-scheme splitting and coefficient functions for electromagnetic deep inelastic scat- 
tering with small x resummations, as given by the NLL BFKL equation with running coupling and 
approximate NLL resummed impact factors. The resummations are combined with a NLO fixed order 
expansion, and the improved quantities thus obtained are stable at small x and significantly suppressed 
with respect to LL results. These results are implemented in a global fit to DIS and related data, and 
the results compared to a NLO fixed order DIS-scheme fit and with a previous LL resummed fit. The 
NLL resummed fit quality is excellent, and constitutes a marked improvement over the purely fixed order 
approach. The input gluon, at Q 2 , = 1 GeV 2 , obtained from the resummed fit is positive and slightly 
increasing as x — > 0, in contrast to the result obtained at fixed order. A resummed prediction for the 
longitudinal structure function Fl is presented. It is positive definite and growing with energy at low x 
and Q 2 , where the fixed order results show a significant perturbative instability. 



1 Introduction 

Current and forthcoming particle collider experiments involve very high energies, such that the momentum 
fractions x of initial state partons are extremely small. For example, the HERA data on the proton structure 
function F2 extends to x > 2 x 10~ 5 QQH]. At these values of the Bjorken variable, the DGLAP splitting 
functions used to evolve the parton distributions, and the coefficient functions used to relate the partons to 
measurable structure functions, are potentially unstable due to large logarithms of the form x _1 ag log m (l/x) 
with n > m + 1 which threaten to undermine the perturbation expansion ordered in fixed powers of as- 
In principle one is able to resum these terms by using the BFKL equation pj], an integral equation for the 
unintegrated gluon 4-point function whose kernel is currently known to next-to-leading order [4l [5] . 



There are several sources of evidence that small x resummations may be necessary when comparing QCD 
with current scattering data, aside from instability in the coefficient and splitting functions for F2 and Fl (the 
fixed order results up to NNLO, and to NNNLO for the F2 coefficient functions, can be found in [5J). Firstly, 
NNLO global fits seem to benefit from the addition of phenomenological higher order terms involving powers 
of log(l/x), whose coefficients are determined by the data [7|. Secondly, the longitudinal structure function 
obtained from the reduced cross-section measured at HERA appears to be inconsistent with the theoretical 
prediction using NLO QCD at small x [8 ], indicating the importance of higher order contributions. We 
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recently carried out a global fit to scattering data at LO in the QCD expansion, supplemented by LL resum- 
mations inclusive of running coupling corrections [9] following the approach of [10] . The BFKL effects were 
seen to significantly improve the description of the low x data when compared to a standard NLO MS-scheme 
global fit. However, the description of the data at moderately high x was not good due to a strong prevalence 
of the resummation effects in this region where the DGLAP theory with no modification ought to be reliable. 

The purpose of this paper is to extend the approach of [10l [9] to NLL order in the resummation, and to 
implement the resummed splitting and coefficient functions alongside a NLO QCD expansion. We will see 
that NLL effects from the BFKL kernel and impact factors suppress the small x divergence in the resummed 
results. This, together with the correct high x behaviour from the NLO expansion, leads to an excellent 
description of the data. Whilst the NLL BFKL kernel is known exactlj{f|, the NLL resummed impact factors 
coupling the BFKL gluon to the virtual photon are not (see [Til HH E2 HH El US] for work in progress). 
These are in principle needed to define the NLL results for the splitting functions P qg , P qq as well as lon- 
gitudinal and heavy flavour coefficient functions^. However, the LL impact factors with the imposition of 
the correct kinematical behaviour of the gluon were calculated in [18j . and in [19j were shown to provide a 
very good approximation to the true NLL impact factors using a comparison with known results from the 
fixed order expansion. In [5U] the exact kinematics calculation was extended as far as possible to the case 
where the virtual photon couples to a heavy quark pair, thus providing all the ingredients necessary for an 
approximate NLL analysis of scattering data. 

The paper is laid out as follows. In section 2 we recall the method of [10] for obtaining small x resummed 
splitting and coefficient functions with running coupling corrections at NLL order. The method is modified 
somewhat from that paper, and so the discussion in this paper is intended to be self-contained. We compare 
our results where possible with alternative approaches [HIES], an d m section 3 we consider the generalisation 
of the variable flavour number scheme of [9] for dealing with heavy flavours to NLL order. In section 4 we 
discuss the details of the global fit, and compare the results of a NLL resummed fit with a NLO DIS-scheme 
fixed order fit. In section 5 we discuss the resummed prediction for the longitudinal structure function, and 
finally our conclusions are presented in section 6. 



2 The BFKL equation at NLL order 
2.1 Fixed coupling solution 

The NLL BFKL kernel is presented in [1] together with the solution of the BFKL equation with fixed as- We 
briefly review this here in order to introduce our notation and also facilitate the comparison with the running 
coupling case to be discussed shortly. First one can introduce the unintegrated gluon density f(x, k 2 ), related 
to the DGLAP gluon g(x, Q 2 ) by: 

/<-«■> -Sfc? 2 - « 

It is convenient to work in Mellin space with respect to the Bjorken x variable, so that convolutions in x are 
unravelled to form products. We use the definition: 

M N [f(x,Q 2 )]=f(N,Q 2 )= [ dxx N f(x,Q 2 ), (2) 

Jo 

3 Contributions to the kernel arising from heavy quarks are yet to be calculated. 

4 We are using the convention of [9] as to what constitutes NLL order, as opposed to e.g. |17| . In our approach, LL order 
for a particular quantity consists of those terms without which there would be no resummation e.g. a^ +1 log n— 1 (l/x) in 
xP qg (as, See [9] for a discussion of this point. 
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where rather than add a tilde to denote the Mellin space quantity, we denote the arguments of each function 
explicitly. Then the BFKL equation can be written schematically as: 



Nf(k 2 ,Q 2 )=Nf I (Q 2 )+a s (k 2 ) / dk' 2 K (k 2 ,k' 2 ,Q 2 ) + as(ti?)}Ci(k 2 ,k' 2 ,QZ) f(k' 2 ), 



(3) 



where as = 3as/7r and the x dependence of the gluon is implicit. The function /i(Qq) is the non-perturbative 
initial condition, and we have taken the gluon at the bottom of the BFKL ladder to be off-shell by an amount 
Qq. The quantities {K, n } are the coefficients of ctg in the expansion of the BFKL kernel. There would be 
an additional term cx a| ln(fc 2 /^ 2 )/Co on the right hand side, but this can be omitted by choosing ^ 2 = k 2 
in the overall power of as- One may deal with the convolution in transverse momentum by introducing the 
further Mellin transform: 

/(iV, 7 )= / dk 2 (k 2 )- x -y{N,k 2 ), 



(4) 

adopting the conventional definition for the Mellin variable. By making the (unwarranted) assumption that 
the coupling as in equation can be held fixed, the result for the double Mellin transformed BFKL 
equation is: 

Nf(N^) = Nf I (N,Q 2 )+aslxoh)+asXi(j)]f(N,l), (5) 

where Xi(l) is the Mellin transform of the kernel ICi(k 2 , k 2 ). Thus the fixed coupling BFKL equation reduces 
in double Mellin space to an algebraic equation, which is easily solved to give: 



/(AT, 7) 



Nfi(N, Ql) 



N -aslxoh) +asXi(l)}' 
The behaviour in a;-space is given after the inverse Mellin transformation: 



f(x, 7) 



1 

2m 



dN 



Nx- N f!(N,Q 2 ) 
N-asxil) 



(6) 



(7) 



where C is a vertical contour to the right of all singularities in the complex N plane, and x(l) is the 
BFKL kernel truncated at the required order. Assuming the perturbative pole given by the denominator of 
equation (|7|) dominates over possible non-perturbative singularities in the initial condition fi(N, Qq), one 
has the leading behaviour as x — > 0: 



f(x, 7 ) = N fj(N ,Q 2 )x- as ^ 



with iVo = asx(l)- The behaviour in momentum space is obtained after a second inverse Mellin transfor- 
mation, which may be written as: 



f(x,k 2 ) = 7T~ / djexp 
2ni J c , 



1 k 2 
asxil) In - +7 In -^2 



iVn 



(9) 



where C is a contour to the right of all singularities in the 7 plane. For very small x, the first term in the 
exponent will dominate, and in the LL case (x(7) = Xo(7)) one may evaluate the integral in equation © 
using a saddle point approximation to obtain the well-known result: 



f(x, k 2 ) - x 



-ias In 2 



k 2 



exp 



log 2 (fc 2 /^) 
56a s C(3)ln(l/a;) 



(10) 



Thus, modulo logarithms in the normalisation, the BFKL gluon density has a power-like growth in x at LL 
order where the exponent arises from the minimum of the kernel Xo a t 7 = 1/2. Carrying out the same 
calculation at NLL order, the power in equation (TIU|) receives a large correction of opposing sign. For nf = 
(no active quark flavours): 

X (l/2) = Aa s ln(2)[l - 6.47a s + ...], (11) 
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where the ellipsis denotes higher order contributions from the kernel, which will certainly not be known 
exactly in the near future. For phcnomcnologically reasonable values of the coupling, the sign of the power- 
like growth in x of the gluon is changed and hence it was initially concluded that the BFKL expansion is 
unstable. Instead, the saddle point calculation is not reliable. The shape of the LL+NLL kernel x(l) is 
completely different to the LL result, with two saddle points at complex values of 7, rather than a single 
saddle point at 7 = 1/2 f2"3"] . 

There are two main solutions to this problem. Firstly, the large correction induced by the NLL kernel can 
be attributed to poles in 7, 1 — 7 which correspond to collinear logarithms of type a r g log™ Q 2 /Qq (n > m) 
in the splitting functions obtained from the DGLAP gluon density in the limits Q 2 ^> Qq, Q\ 3> Q 2 . It is 
possible to resum these contributions [H] and thus obtain a kernel with a stable minimum at NLL order 
at 7 = 1/2. Secondly, the calculation above ignores the fact that the QCD coupling as runs with energy, 
which must be taken into account beyond LL order in the resummation. This changes the nature of the 
BFKL equation so that it is a differential equation as opposed to a purely algebraic one. We will sec in this 
paper (as has already been observed in |10j ) that in solving this equation one is no longer concerned with 
the behaviour of the kernel at 7 — 1/2. Alternative approaches to small x resummation of splitting functions 
|21|.[22j involve combining the resummation of collinear poles with inclusion of the running coupling. We do 
not find it necessary to resum the kernel, and consider the solution of the running coupling BFKL equation at 
NLL order in the following subsection. This is precisely because the deep inelastic scenario where Q 2 ^> Q\ 
focuses on the region 7 = 0. Single scale processes, where there is a large scale at both ends of the gluon 
ladder, are a different matter. 



(12) 



2.2 Running Coupling Solution 

Taking ag(fc 2 ) as the coupling in equation ©, one may substitute the LO expression for the coupling: 

a s (k 2 ) _ 1 
4tt /3 ln(fc 2 /A 2 ) 

(where (3q = (11 — 2n//3)) and multiply through by ln 2 (fc 2 /A 2 ) before Mellin transforming in transverse 
momentum to obtain the second order differential equation: 

d 2 fh,N) d 2 fj(j,Q 2 ) 1 rf(xo( 7 )/(7,AQ) , tt v ,,., m n ^ 

introducing /3q = 7r/3o/3. A comment is in order regarding the use of the LO running coupling. Given that 
our aim is to build a NLL resummation upon a NLO fixed order expansion, one should in principle use the 
NLO running coupling in the BFKL equation. This leads to a more complicated equation in momentum 
space, which is not diagonalised by a Mellin transform. Indeed, it is not clear how an analytic solution can 
proceed in this case. The use of the LO coupling can be justified given that differences between the LO and 
NLO couplings can be absorbed by a change in A. We discuss in more detail later how to combine resummed 
quantities with the fixed order expansion. 



Equation (fT3|) may be solved by making the ansatz [25] : 
where: 

H 



f(N, 7) = exp Ml) ex P (^M ) dry, ( 1 4 1 



Xt = J i X NLo(l,N)dj, (15) 

and xnlo and ^(7) are to be determined. Substituting this into equation (| 13[) and choosing ^(7) via: 

dA(j) _ A(j)\xnlo(% N) - X o(7)] rf 2 //(7) 



g?7 f3 N d-y 2 



(16) 
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one finds: 



PoN dx ^° - {xnlo ~ Xo)xnlo = N[f3 x'o(j) ~ Xi(7)L 



(17) 



where a factor of ir/3 has been absorbed into Xi(l)- This is a first order nonlinear differential equation for 
Xnlo- It may be solved iteratively by expanding Xi(l) m powers of TV (the u uo expansion" of |25j in their 
notation): 



One then gets: 



Co (7) 

a (7) 
6(7) 



Xnlo 



xo(t); 

Xi(t). 
Xo(7)' 
1 

Xo(7) 



(7,^) = E^(7)^ J 



(18) 



/ Xi(7) 
VXo(7) 



Xi(7) 
Xo(7) 



(19) 



Note the term involving (3$ in £2- This is the first of a series of terms 
in the iterative solution of equation (|17p from: 



Pod, 
Co ^7 



/3g 1 A™ in the TV-expansion, arising 



(20) 



where the ellipsis denotes sub-leading running coupling corrections. A decision must now be made as to 
what constitutes NLL order. It could be argued that one must keep the whole of equation (fT5)) in defining 
the NLL gluon density, which arises from a complete solution of the NLL BFKL equation. However, this is 
impossible to achieve in practice due to an infinite number of terms. Instead, one might consider truncating 
at 0(A) and keeping the set of leading running coupling corrections given by equation (|2"0)) . Again one has 
the problem of infinitely many terms, and in fact they should not be included, which can be seen as follows 
[25] . Expanding in 7 each of the coefficients of the A-expansion one has: 



n=l 



7 



and we will see that this leads to a gluon density of form: 

with t — log(fc 2 /A 2 ). This will lead to a contribution to the LO anomalous dimension P gg of form: 



p(0) 
99 



a s 1 



KN 7 



(21) 



(22) 



(23) 



where if the above-mentioned terms are kept, k n oc /3g _1 . We know, however, that the LO anomalous di- 
mension does not contain higher order terms in /3o, and so the terms of form /3q _1 A" in the A-expansion 
must somehow be cancelled. To see how this works one must consider the BFKL equation at higher orders 
in the high energy expansion. If one expands the kernel to NNLL order, the BFKL equation will become a 
third order differential equation in Mellin space due to the extra power of as in the kernel (assuming the 
LO coupling is still used). The Mellin transformed NNLL kernel X2(l) contributes to £2(7) in such a way 
as to cancel the term in (3$. Similarly, the Mellin space BFKL kernel at N"LL order in double Mellin space 
is an (n + l) th order differential equation. The higher order kernels contribute to the coefficients of the 
A-expansion in such a way as to cancel the running coupling terms discussed above, as must be the case. 
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Thus, we now have the gluon density given by equation (|14[) . where the series for xnlo is truncated at 
O(N). For I7I < 1 one may expand X\ as: 



/ 00 

Xi( 7 ,A0 = (\og 1 + lE +^a„7 2 ™ +1 ) -iV Qlog 7 + c + ^ 



Cn7 ' 



(24) 



where 75 is Euler's constant, a n = 2((2n + l)/(2n + 1) and the {cj} follow from the NLL BFKL kernel 
[10] . Thus the integral on the right-hand-side of equation ([14"]) has behaviour ~ 7 



(l-c,ftJV)/ftJV 



Hence, an integral of this from to 7 behaves like ~ 7 



l + (l-c,ftJY)/ftJV 



as 7 



0. 



cancelling the singularity of the 



prefactor. The leading singularity of the gluon from this contribution is then at 7 = —1, and so up to 
power-suppressed corrections ~ A 2 /Q 2 in momentum space one may shift the lower limit of the integral 
from 7 — > leading to factorisation of the gluon density into a perturbatively calculable piece (the prefactor 
in equation (|14p ) and a non-perturbative piece which ultimately is combined with the bare gluon density in 
the proton. Transforming back to momentum space, the perturbative piece of the integrated gluon density 
is: 

^ /-1/2+10O r/3 

Oh(N,t) = — ^exp[ 1 t-X 1 ( Jl N)/(p N)]d 1 (25) 

where the inverse power of 7 comes from considering the integrated gluon, and we have used the same 
notation as [10]. The factor f"° in equation (|25|) is given by: 



fftj _ 



exp 



1/2 



(26) 



and corresponds to the contribution in the exponent of equation (|25[) from the running coupling contribution 
to the NLL kernel, effectively providing the corrections required from making the simple choice of scale k 2 
in the coupling. We factor this out explicitly following [10]. From now on, the notation Xi(t) refers to the 
kernel without this contribution. 



2.3 Solution for the BFKL Anomalous Dimension 



The solution for the perturbative piece of the integrated gluon density of equation (]25[) can be evaluated 
by numerical integration. Alternatively, one may obtain an analytic expansion in momentum space by 
substituting the expansion of equation (|24| into equation Q25p and deforming the Mellin inversion contour 
to enclose the negative real axis. One obtains (taking into account the discontinuity across the cuts on the 
real axis): 



G E (N,t) =-sin 



x exp 



n(l-ciN) 



exp 



It- 



Substituting y — 7 t and using the result: 



a„7 



1e + c \ r° 
2n+1 +Nc nl N ) 



j/3 ( 7 ) 7 -(l-c 1 A')/(ft>JV)-l 



dyy x 1 exp( 2 /) 2 /" 



(-l) A - 1+, T(A + n) 
(-l) A - 1 A„(A)r(A) 



(27) 

(28) 
(29) 
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after Taylor expansion of the integrand, one obtains: 



Ge(N, t) = sin (« l r«W) exp ( -12^*») T ( Jlz^l) & 

f OO 

x i+E 



k n=l 

OO OC 



i„(l/(/3 iV))4 C " (1//3o) 



r"A. 



l-c,A T 



with the coefficients A n and C rn dehned by: 



f + £ A n (l /(P N)h n = f° exp (Jjj J2 a ^ n ) 5 

oo / 1 oo \ 

f + £ ^UV^oh" 1 = exp U- £ c„ 7 " , 

m=l \"° n=l / 



(31) 
(32) 



The i-independent pre-factors in equation (|30[) can be absorbed into the non-perturbative gluon density. 
From equation (j2T)|) one may write for the functions {A„} defined in equation (f2"5)l : 

n-l 

A„[(l - C/ 7V)/(/3 ^)] - ]T (-l) m ^-(/5 ^V)- n+m , (33) 

with do, n = 1 and <i m , n +i = ^m,n + n d m -i,n- Setting m = n in this recurrence relation yields d n .n+i = 
ndn-x,n and hence: 

d„_i,„ = (n-l)! ! (34) 

and so the series in equation (f3T)]) has factorially divergent coefficients at high n. Thus the series is asymptotic, 
and must be truncated^] at some finite order n = uq. In the LL case of |f 0j . no = 5 was found to be sufficient. 
There are other choices that must be made regarding what is meant by NLL order, and in [TU] it is noted 
that equation (|30|) contains NNLL terms coming from the following origins: 



• The C m (Po) can be expanded as a power series in (3 1 where terms of O(/3 2 ) ar e really NNLL 
contributions. 

• The A„((l — ciN) / ((3qN)) functions in equation (|30|) can be expanded as a power series in TV, where 
terms 0(N 2 ) are NNLL contributions. Using Taylor's theorem one may write: 

( (l-c t N) \ _ (_ 1 \ c, dA n (-l/p N) 

where the ellipsis denotes NNLL terms. In [10] these were discarded, although this is not strictly necessary 
as one is free to modify a NLL resummation up to NNLL terms. Indeed, making the truncations advocated 
in [TUJ leads to resummed splitting functions that are significantly different from the fixed order NLO results 
at moderate x, particularly P qg which shows a pronounced dip 26J. Indeed, consideration of the anomalous 
dimension j + shows that this truncation is incorrect, as described below. Another advantage of keeping the 
full arguments of the A„ functions is that this argument is very small at N = 1. When nf = 4 for example: 

-^T" = -^(1- 0.9417V) = 0.08, N = l. (36) 



3 The series is, however, summable in principle due to the oscillatory factor (— l) m in equation )33t . 
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The splitting functions obtained from the gluon density will also have a small first moment, and hence 
satisfy the momentum sum rule to a very good approximation. Thus we choose to define the gluon density 
by equation (|30|) with no further truncation. The BFKL anomalous dimension (inclusive of running coupling 
effects) is now given in principle by: 

1+ - NLL ~ GUN,t) at ' (37) 

which corresponds to the leading eigenvalue of the singlet anomalous dimension matrix: 

NLL „, , r „„s 
7+ = Igg + ~cpi9 V™) 

at high energjU. There is now a further choice that has to be made in defining the anomalous dimension. 
Decomposing the gluon into LL and NLL part^ as Q\ = Q^ L + G^ LL , one has: 

1 <5 \ r LL , r NLL] 

1+ - NLL = W 7 +G¥ Z di [ E E 1 

C NLL\- 1 i a 

^%r) s|t| + (39) 

The question now arises: should one keep the whole of the first factor on the right-hand-side, or expand it 
to NLL order? The latter choice yields: 

■^-(l-^gklpi' + tfH, (-0) 

and is a further truncation in addition to those already described in the gluon density. The poles in equation 
(jl(7)l now only come from the zeros in Qj^ L and this ultimately leads to a corrected splitting function behaving 
like: 

P+ lW ii~a:''[l-A7 7 log(l/a;)] (41) 

rather a shifted power-like behaviour: 

P+.NLL^X^, (42) 

as argued in [lOj . The latter behaviour is physically more sensible, and the argument can be made more 
compelling by examining the splitting functions obtained using equations (|39l40p . These are shown in figure 
[TJ where the truncated splitting function in this plot corresponds to a full removal of NNLL and higher terms 
(apart from running coupling corrections) after using equation (|40[) to expand the gluon denominator. There 
is a marked difference in asymptotic behaviour. The NLL result with no truncations dips gently under the 
NLO fixed order result, ultimately rising at low x. The truncated expression decreases far below the fixed 
order result and never rises at small x, consistent with a shift in behaviour given by equation (|41[) . Thus to 
obtain the correct physical behaviour for the structure function, one must use the complete denominator of 
equation (|3"T|) . 

The intermediate possibility exists of truncating the solution (|30|) for the gluon, but still using the complete 
gluon density in the denominator for the anomalous dimension. However, this introduces spurious terms 
involving inverse powers of /3o, even at orders below that corresponding to the truncation of the gluon density 
in 7-space. Such terms are non-physical and can be avoided as far as possible by not truncating equation 
(|30| . The Pq 1 terms introduced by expanding the denominator ultimately cancel with those obtained by 
keeping the full (3q dependence in the C' m (Po) coefficients discussed previously. 



°The alternative approaches |21U22j present results for 7 39 . However, these are in the limit rif — > where one has 7+ — > 7 99 . 
7 The NLL piece contains some higher order terms than strictly NLL order as discussed after equation (1306 - 
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Figure 1: Power-series result for the BFKL anomalous dimension P+^nll for t = 10, n/ = 4 obtained from 
using the full gluon solution of equation (pJO")) . together the result obtained after truncating to formal NLL 
order and the asymptotic NLO (fixed order in as) result. 
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There is a power-suppressed correction to the anomalous dimension obtained from equations (|30I37[) . arising 
from the fact that the expansion of Xi("f, N) is only valid for 7 > — 1 and also from the truncation of the 
asymptotic series for Xi^Jf]. To evaluate this correction, one must evaluate the anomalous dimension from 
numerical integration of equation (|25|) . For this purpose the gluon density can be rewritten as: 



1 />l/2 + 200 

Gh(N, t) — — / djf° ( 7 ) 7 - 1+c '/^ exp 

Jl/2-joo 

n 



7* -^o(7) + ^i(7) 



(43) 



where A (7) = / 1 / 2 Xo(7)d'7 and ^1(7) is that part of ^1(7) that gives J2 n =i c "7™ when Taylor expanded. 

A problem arises in implementing this correction, in that a closed result for Xi("f) does not exist as one 
cannot integrate £1 (7) of equation (|19p analytically. Given that the correction is found numerically, however, 
it is sufficient to parameterise £1 along the integration contour. We use the contour: 

f0.4-(i-lV, -oo<u;<0 . , 

1 \ 0.4+0-l)w, Q<w<oo ' K ' 

which is a deformation of the Mellin inversion contour that encounters no singularities. In practice it is 
sufficient to parameterise £1(7) for < w < 3, and then one has along the upper branch of the contour: 



P w 

X l {w) = {I-l) / &[0.4+(i-l)t 
Jo 



]dw. (45) 



Then Xi(w) is given along this contour by: 

Xi(w) = Xx{w) - c ; [log(0.4 + (t - l)w) - log(0.4)] - c , (46) 

where: 

co = / 6(7)^7 (47) 

JOA 

which can be integrated numerically. After evaluating Q\ and dQ Y E jdt, one generates the anomalous dimen- 
sion according to equation (|37p . We find that the correction to the power series result is very large. It can 
be stabilised somewhat by taking fewer orders of 7 in the expansion of the anomalous dimension in double 
Mellin space (i.e. no — 4 instead of no = 5), but is still comparable with the power-series result for t < 11. 
Thus, instead of parameterising the correction as in |10j . we choose to parameterise the anomalous dimension 
directly at lower t values. Using the integral representation of equation (|43|) . one can obtain values of the 
splitting functions for particular values of N and t and fit these to a function of form: 



±77 



J- '"U j 

P(N,t)= Ep^^T, (48) 

m— — 4 n=0 

where negative powers of t happen to achieve a better fit for P+ and positive powers improve the fit for other 
quantities to be discussed in the next section. Several ranges of t are needed to achieve a good fit, and one 
must also calculate separate parameterisations for each value of n/, due to the dependence on the number 
of flavours in the NLL BFKL kernel. At high enough t one may use the power series result for the splitting 
functions. The power series results must be expanded to many orders of as (around 20) and due to their 
length are not presented here. 



2.4 Inclusion of Impact Factors 

In order to complete the phenomenology at NLL order for massless quarks, one must calculate the resummed 
splitting function P qg . One also needs the coefficient function Cl s to obtain a prediction of the longitudinal 

8 This is also the case at LL order with running coupling. See 1101 . 
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structure function. Thus one needs the impact factors coupling the virtual photon to the gluon evaluated 
at NLL order. These are not yet available but one may proceed by using the LL factors calculated with 
the imposition of exact gluon kinematics [18j . shown in [19j to be a good approximation to the complete 
higher-order results. Their interpretation in terms of physical quantities stems from the /cy-factorisation 
theorem [27J[S5], which in double Mellin space dictates for the longitudinal structure function: 



F L (~,,N) = h L (j,N)g(j,N), 



(49) 



where the gluon density 5(7, N) factorises further into perturbative and non-perturbative contributions as 
discussed above. From equation (|49p one may straight-forwardly associate the impact factor hh with the 
longitudinal coefficient function Cl iS in double Mellin space, where this coefficient is in the "Qo" factorisation 
scheme where the gluon is obtained by solution of the BFKL equation with an off-shell non-perturbative 
initial condition^. The case of F 2 is not so straightforward. The quantity analogous to hh in equation (|49|) is 
the sum of impact factors /it(7, N)+hi J ('y, N) arising from transversely and longitudinally polarised photons 
respectively, and this has the leading behaviour 7 -1 as 7 — > 0, corresponding to a collinear divergence of 
the momentum space impact factor. Put another way, the structure function F 2 is proportional to the (non- 
perturbative) quark singlet distribution at LO in the QCD expansion, and one does not expect perturbation 
theory to be able to describe this behaviour. Instead one may consider the quantity: 



dF 2 



dlogQ 2 



h 2 (j,N)g(j,N), 



(50) 



which starts at 0(as) in perturbation theory. The left-hand side denotes the double Mellin transform, 
and equation (|50[) serves to define the impact factor h 2 . The interpretation of this quantity depends on 
the factorisation scheme, and in general represents a mixture of the coefficient function C 2 . g (j,N) and the 
anomalous dimension 7 9g . If one chooses the DIS factorisation scheme [35] , one has the simple identification 
of h 2 with 7^ /s in double Mellin space. 

The exact kinematics impact factors of |18j can be expanded in N: 

h a ( 1 ,N) = h^( 1 ) + Nhi 1 \ 1 ) + 0(N 2 ), 
corresponding to truncation at NLL order. Then equation (|43|) can be generalised to give: 



(51) 



1 /•1/2-Hoo 

Fh(N,t) = — / d 7 M7,A0/ & (7)7- 1+Ci/ ^exp 



2m 



l/2-ioo 



(52) 



Proceeding analytically as for the gluon density one obtains for the structure function: 

00 

Fh(N,t) oc *(i-c.iv)/(/3oiV) J- ^(o) + Nh (i)- 



r=0 



00 00 



A n (l/(f3 N)) 



0o 



EE 



A n (l/(PoN))C m (l//3o) 



n—1 m—1 

where the impact factors have been expanded as follows: 




00 

/ la ( 7 ,^V) = ^[e ) +iV/ 1 W 



7 • 



(53) 



(54) 



9 The "Qo scheme" is widely used in the literature to denote a gluon arising from solution of the BFKL equation, whether or 
not an off-shell regularisation is used in practice. It is thus a type of scheme rather than a specific scheme, due to the number 
of choices to be made when solving the BFKL equation. 
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Note that there are NNLL contributions in equation (|53[) arising from the product of the O(N) piece of 
the impact factor with NLL terms from the gluon. This is not a problem in that they are beyond the 
NLL nature of the resummation. Truncation of the structure function results, as in the purely gluonic case, 
in unphysical terms involving inverse powers of (3q. In order to avoid such contributions, I choose not to 
truncate equation (|53p . Once again the full argument of the A„ functions is kept, which will ensure that 
momentum is conserved by the resummed splitting functions to a very good approximation. 

Note that the inclusion of impact factors has led to terms ~ 0(N) in the structure functions, which conse- 
quently lead to similar terms in the quantities P qg , Cl 3 in addition to 0(N ) terms which are already present 
at LL order (see QUI)- In x-space, these correspond to the Dirac function 8(1 — x) and its first derivative. 
However, implementation of #'(1 — x) is numerically cumbersome and is somewhat unnatural given that 
there are no real terms proportional to <5'(1 — x) in the fixed order splitting and coefficient functions. The 
0(N) terms are an artifact of expanding about N = in the resummation. Instead of implementing them 
as Dirac functions, one may choose a suitable function whose expansion about N = has similar properties. 

In more detail, consider a resummed splitting or coefficient function that has form: 

P = 5 1 (t)N + S (t) + O(N- 1 ). (55) 

One wishes to replace the Dirac terms of P with a function whose Taylor expansion about N = reproduces 
these first two terms. Furthermore, it is important to preserve the first moment (value at N=l) of a splitting 
function to ensure compliance with the momentum sum rule. We have already seen that the power series 
results for the splitting functions ought to have a very low first moment, and this property must be preserved 
when implementing the Dirac terms. One thus has three pieces of information with which to construct a 
function. There is relatively little sensitivity to the choice, and a suitable form which has no 0(N~ 1 ) 
component ij^l: 

P'(x) = ^P- + k 2 (t) + k a (t)V^, (56) 

y/X 

which in TV-space becomes: 

fci(t) k 2 (t) k 3 (t) 



P'{N) = 



AT + 1/2 N + l N + 3/2 



= f 2fci(t) + k 2 (t) + Hfc 3 (t)) + U4fci(t) - k 2 (t) - ~fc 3 (t)J N + 0{N 2 ). (57) 

The need to reproduce the first two terms of P, together with the first moment condition then gives the 
simultaneous equations: 

2fci(t) + fc 2 (t) + |fc 3 (t) =<*(,(*); 
-4fci(t) - k 2 (t) - ^(t) = «i(t); 

|*i (t) + \k 2 {t) + h 3 (t) = S (t) + 8i(t), (58) 
3 2b 



which one may solve to give: 



fc2(t) = -8<f (t)-24ai(t); 

k 3 (t) = ^-[26 (t)+5S 1 (t)}. (59) 



'An alternative form e.g. P'(x) = fci(t) + &2(t) log(l/a;) + kg(t) log(l/a;) 2 gives similar results. 
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If one simply omits the 0(A) terms in the resummed coefficient and splitting functions, it is not possible to 
achieve a good fit to scattering data as the resulting splitting functions lead to a large violation of momentum 
conservation in the evolution. 



One must again parameterise the splitting and coefficient functions for t < 11. A function of the form 
is sufficient, with positive powers of t for both P qg and Cl s - We have so far ignored those quantities which 
couple to the quark singlet distribution i.e. P gq , P qq and C^q- At LL order, these are related to the analogous 
gluon quantities via the colour-charge relation: 



q ~ C A 



F - p(°) 



(60) 



where F a is a LL resummed coefficient or splitting function coupling to singlet species a, and its leading 
order contribution to the fixed order expansion. We also use this relation to define resummed quark singlet 
quantities at NLL order. This is not formally true at this order, but is a reasonable approximation with a 
similar uncertainty to that arising from using the exact kinematics estimates of the NLL impact factors. 

Before presenting results for these quantities, the resummation from the BFKL equation must be combined 
with the NLO fixed order expansion. This poses a potential problem in that the gluon in the Qo-scheme is 
in principle different to that in the fixed order DIS scheme. We discuss this in the next subsection. 

2.5 Transformation from the Q Scheme 

The resummed splitting functions found so far are in the Qo scheme, where the gluon is defined via solution 
of the BFKL equation. However, the gluon in the DIS scheme at fixed order is given in terms of parton 
densities in the MS scheme, which arise from using dimensional regularisation. Thus the two gluons are not 
the same, and in principle a scheme change is needed before one can add the resummed quantities to their 
fixed order counterparts. 

Let us examine this in more detail by looking at the anomalous dimension given by equation (|37[) . which 
when expanded out in powers of as in A-space gives: 

7+ JVLi = a s (Jj - 0.91667 - 0.0061728n/J + af ( - °- 10 ^ 7n / + 0.09713671/ 

+ 0.00065412rr^ + 0(N). (61) 

One may compare this with the relevant fixed order quantity, which is 7 99 + (4/9) 7gff as discussed previously. 
The NLO DIS scheme result, expanded about N = 0, gives: 

7+ =a s - .91667 - 0.0061728n/j +a|(- • 105 ^ 7n f _ .H579n/ + 2.1864 

+ 0.0058299n^ +0(N). (62) 

The resummed and fixed order results are not equal, and this is partly due to the difference in gluons in 
the resummed and fixed order framework^]. The solution is to implement a scheme change in the kx 
factorisation formula, making the ansatz in double Mellin space [30] : 

S( 7 ,A) =i?( 7 ,A)£ DIS ( 7 ,A), (63) 

where G Dls is the fixed order unintegrated DIS scheme gluon density and R(j, N) is a perturbatively calcu- 
lable transformation factor. A fully consistent NLL resummation consists of implementing the BFKL kernel 



11 There are also NNLL contributions, distinct from the running coupling corrections, not predicted by the resummation. 
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as above, but also including an all orders in 7 result for R(j,N) truncated at 0(N). The O(N ) part of 
i?(7,iV) is known [31], but the O(N) piece is not. It would be found by formulating the BFKL equation 
using dimensional regularisation, and solving with running coupling. This does not seem to be analytically 
tractable at present [32], and thus a compromise must be reached. One might instead try exploiting the fact 
that differences between the fixed order and resummed gluons can be neglected at orders of as beyond the 
extent of the fixed order expansion, and thus expand the i?(7, N) factor as follows: 

JZ(7, N) = r + ri7 + 0{j 2 ) + [r 2 + r 37 + 0{j 2 )]N + 0{N 2 ). (64) 

This corresponds to truncating the NLL resummed scheme transformation factor at NLO in the fixed order 
expansion. The generalisation of equation (|25[) to include the R factor is straightforward, and in practice 
it mimics the addition of an impact factor. One may find the coefficients Ti by comparing the anomalous 
dimension jbfkl obtained by inserting equation into equation (j2"5)) with the fixed order 7 + . However, 
it is easy to see that ro = 1 and r 2 = from the fact that at LO the scheme transformation must be trivial 
(R = 1). Also n = as the LL factor #(7) = 1 + C(7 3 ) [31 ]. Thus up to NLL and NLO orders, one has 
i?(7, N) — r^jN and the resulting anomalous dimension turns out to be: 

7+JVLL =a s (^- 0.91667 - 0.0061728n/j + a 2 s ( - °- 10 ^ 7n / + 0.097136?!/ 

+ 0.00065412?^ - /3 r 3 ) + O(N). (65) 

This can be compared with equation (|62[) , and in principle the coefficient r 3 can be read off. However, 
the term at O(N ) at NLO in equation (JB5J) is a sum of a running coupling correction (predicted here by 
the resummation) and a NNLL term (not completely predicted by the resummation) . It is not possible 
to disentangle the former contribution from the latter, and hence unambiguously determine r 3 . However, 
it seems that one can neglect this contribution to the scheme change to a very good approximation. In 
[19j we have compared the fixed coupling results for the NLL resummed P qg (using the exact kinematics 
approximation) with the exact NLO, and NNLO splitting functions. The O(N ) term in that splitting 
function at NLO and the 0(-/V _1 ) term at NNLO were estimated extremely well - suggesting that the 
scheme change from the Qo scheme to a fixed order type DIS scheme (where partons are dimensionally 
regularised) is a very small effect. It is certainly a less significant effect than using the exact kinematics 
approximation for the NLL impact factors. The issue of transforming from the Qo scheme to the MS scheme, 
however, is still unresolved. 



2.6 Matching to the NLO expansion 

So far we have obtained NLL resummed splitting and coefficient functions, using the BFKL equation and 
the exact kinematics impact factors. In order to use these in a fit to scattering data, one must combine 
them with the NLO fixed order QCD expansion. Care must be taken to avoid double-counting of terms. 
Schematically one has a total splitting or coefficient function given by: 



ptot. _ pNLL , pNLO _ pNLL(0) _|_ pNLL(l 



(66) 



where the bracketed term on the right-hand side subtracts the double-counted terms, and p NLL i n ) jg the 
contribution of the resummed quantity to the n th order term in the fixed order expansion. This is complicated 
in practice by the fact that the resummed quantities are not known accurately as a power-series in ag, and 
must be parameterised according to equation (|4"8)) . Nevertheless, one can use the LO and NLO terms of the 
power series solution (from equations (|53l30p ) for each quantity to determine these terms. As an example, 
consider j qg at t = 6 and n/ =4. This can be found by numerical integration of equations (|27l [52| and 
parameterised by: 

ll+nll ctsnt ( crniAr , , cc , 0-2294 0.03443 0.01381 ^ 
7„™ m . = [ — .5501JV + .2466 + — j^- + ^ \ . (67) 
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Figure 2: The total NLL resummed splitting functions P qg and P+ for n/ = 4 and t — 6, compared with the 
NLO fixed order results. 



However, this includes some of the high energy limit of the LO and NLO anomalous dimension, which must 
be subtracted to avoid double counting. These double counted terms are obtained from the power series 
solution for ^ qg up to NLO, which is: 



LL+NLL = a Snf 
lqg,pow. 



2 16T 

I - FOSSA' o s ( — 6.484 + 4.1847V 



+ 0(a%). 



(68) 



Subtracting (f6"5| from ([6"?]) gives the purely higher order corrections to be added onto the fixed order expan- 
sion. 



Our results for P qg and P+ are shown for nf — 4 and t — 6 in figure [2] and alongside their LL counterparts, 
both with and without the inclusion of running coupling effects, in figure [3] The inclusion of the running 
coupling at LL order suppresses the resummed splitting functions, as previously noted |10j . The NLL 
corrections from the impact factors and kernel then lead to a further suppression of the small x divergence. 
Both the splitting functions dip below the NLO result, before growing at very small x. One must also bear 
in mind the missing Dirac terms in figure [5J and [3J These carry a positive first moment in the case of P qg , 
but negative for P gg . The NLL corrections delay the onset of the asymptotic behaviour until even lower 
values of x. The same splitting functions at a higher value of t = 8 are shown in figured] 

It is interesting to compare our results for P + with those of the ABF [3T] and CCSS [55] approaches, and a 
comparison is given in figure[5] All three splitting functions show the same qualitative behaviour. There is a 
pronounced dip below the NLO result at moderate x, followed by an eventual asymptotic rise as x — > 0. At 
high x the behaviour is governed by the fixed order DGLAP evolution kernel, as required in each approach. 
Qualitatively, the CCSS function has a slightly deeper dip, and we note that both the CCSS and ABF results 
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Figure 3: The NLL resummed splitting functions P qg and P + for nj = 4 and t = 6 (solid), together with 
the LL result with (dashed) and without (dot-dashed) running coupling corrections. Also shown is the NLO 
fixed order result (dotted). 
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Figure 4: The NLL resummed splitting functions P qg and P + for n/ = 4 and i = 8, compared with the NLO 
fixed order results. 
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Figure 5: Comparison of the splitting function in figure (|3|) with the CCSS and ABF results for as = 0.2 
(t ~ 7.5 for the LO coupling). The CCSS and ABF results have nj = 0. Results taken from [21]. 
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Figure 6: The NLL resummed coefficient Cl s at t = 6 and n/ = 4 (solid) together with the LL result with 
(dashed) and without (dot-dashed) running coupling corrections. Also shown is the NLO result (dotted). 



show a slightly higher rise at asymptotically low x than the splitting function presented in this paper. Part 
of this may be due to the fact that both of the alternative approaches involve the rcsummation of collincarly 
enhanced terms in the kernel, which act to suppress the impact of the NLL contribution. We see from the 
figure that this is not a large effect. However, it must be noted that the CCSS and ABF curves in figure 
[5] have n/ = 0. The effect of a higher quark flavour is increase the splitting function at low x, implying 
that resummation of the kernel may in fact be more significant than appears to be the case in figure El It 
would of course be extremely interesting to see a comparison of the quark splitting functions, as well as an 
application of all three approaches to experimental datj^l. 



In figure [5] we show the resummed longitudinal coefficient Cl„ at t — 6 and rif — 4. Once again, the NLL 
corrections suppress the result at small x when compared with the LL result, beyond the suppression already 
induced by running coupling corrections. The results in figures [5J[5] do not include the effect of the O(N ) 
and 0(N) terms from the resummation, which are implemented according to the prescription of equation 
(|56[) . In the example of "f qg at t = 6 and rif — 4 given above, one finds: 



P' = 

19 



3tt 



.6340 



(A +1/2) 



4.879 



4.213 

(iV + 3/2) 



(69) 



Likewise, for Cl.q at this value of t one has: 



C 



3tt 



.4046 



(N + l/2) 



.2850- 



.1945 



(A + 3/2) 



(70) 



12 A discussion of impact factors in both approaches can be found in |33l I34| . There is a preliminary investigation of P qg in 
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Figure 7: The resummed + NLO results for P qg and Cl 9 at t = 6 and nf — 4 where the Dirac terms have 
been implemented according to equations (|69l l70f , compared with the NLO results. 



One can examine the form of the resummed + NLO quantities in x-space. These are shown in figure [7] 
Comparing this with figures [2] and [51 one sees that the Dirac terms have little visible effect on Cl 9 (mainly 
due to the pronounced high x behaviour in the NLO coefficient) whereas in P qg the dip at moderate x is 
made more significant. However, a direct physical interpretation of figure [7] is not really possible. Other 
functional forms for implementing the Dirac terms in iV-space (where they are constrained) can lead to 
quantitatively different shapes in a;-spacJ^l. but similar results in convolutions. 

We now have the complete set of massless resummed splitting and coefficient functions needed for phe- 
nomcnological investigations. Before a fit can be made, however, one must also consider rcsummation effects 
in the heavy flavour sector. This is the subject of the next section. 

3 Treatment of Heavy Flavours 

In [9] we outlined the definition of the DIS(x) scheme for the consistent implementation of small x resumma- 
tions in the heavy flavour sector alongside a fixed order QCD expansion. We refer the reader to that paper 
for details, including a full discussion regarding consistency of the variable flavour number scheme between 
the fixed order and high energy expansions. Here we extend this scheme to NLL order in the resummation 
as far as is possible. 

Below the matching scale Q 2 — M 2 for a heavy quark of mass M, the heavy flavour contribution to the 

13 The NNLO splitting function P qg does, however, display a narrow but deep dip at high x, as can be seen in 1 191 . 
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structure function F2 is given by the fixed flavour description: 



F H - 

r 2 ~ 



CiI' LL (l, N, M 2 /Q 2 ) + C£f ""(7, N, M 2 /Q 2 ) 9 <«%, N). (71) 



y 2,g Ui"i™ I^C J-r ^2,g 

Comparing with the fc^-factorisation result: 

F* = h 2 ("f,N,Q 2 /M 2 )g(N n ), (72) 

one identifies the fixed flavour coefficient with the impact factor ^2(7, N, Q 2 /M 2 ) coupling a virtual photon 
to a gluon via a heavy quark-antiquark pair. A calculation of this quantity to NLL order will certainly not 
be available in the near future, and by analogy with the massless case one may use the LL factor calculated 
with exact gluon kinematics which has been presented in [20) . 

The factor h diverges as Q 2 /M 2 — > 00 to give the collinear singularity of the massless result, and so when 
Q 2 > M 2 one must use the NLL variable flavour description: 



F 2 H = q H , 



CH' LL h, N, M 2 /Q 2 ) + CY Q NL \1, N, M 2 /Q 2 ) g*f+V, (73) 



J 2,g Ul 1 '! 1 " ) ~r ^ 2 ,g 



,(«/+!) 



where qn+ = Qh + <1h is a parton distribution for the heavy flavour. The (rif + l)-flavour partons are related 
to their n /-flavour counterparts by resummed heavy matrix elements Ay , which at NLL order are contained 
in the relational: 

q H+ ( 7 ,N,M 2 /Q 2 ) = [A^( 7 , A,M 2 /Q 2 ) + A^ L ( 7 , N, M 2 /Q 2 )] g&>\ (74) 
g( n f+ 1 \ 1 ,N 1 M 2 /Q 2 ) =Af g LL ( 7 ,N,M 2 /Q 2 )gM. (75) 

Substituting these in equation (|73[) and equating coefficients of the gluon at LL and NLL orders, one finds: 

n VF,LL _ n FF,LL a LL . / 7r n 
°2,g — °2,g A Hgi \<°) 

r ,VF,NLL _ n FF,NLL aNLL n VF,LL aNLL (~*\ 
°2, S -°2,g A Hg °2,g A gg • {" ) 

Equation ([76]) has already been derived in [9], and equation ([77]) is its NLL equivalent. The equality must 
be true for all values of Q 2 /M 2 . Taking Q 2 /M 2 — ► 00 gives: 

aN LL n FF,NLL\ / 7q n 
A Hg = C 2,g I.SL_oo" (78 ^ 

M 2 

where we have used the fact that C^ g — > as Q 2 /M 2 — > 00 in the DIS(%) scheme. A problem now occurs 

111 that A% g LL and C^ g are both still unknown in equation (|T7[) . One needs a further equation to fully 
define the variable flavour scheme, and this must come from a definition of A gg . Unfortunately, this is not 
possible at present, as in principle A gg is determined from a calculation of the gluon density with heavy quark 
mass effects included. This would involve knowledge of the NLL BFKL kernel with heavy quarks included 
in the fermionic contributions, and this is not available. One can only neglect the term involving A gg and 
hope that that it is not significant!^!. We thus define the (approximate) NLL variable flavour coefficient to 
be: 

n VF,NLL n FF,NLL aNLL 
°2, S - ^2,3 ~ H.g 

r,FF,NLL ^FF,NLL i lnr .\ 
= C 2,g -C2,g loi^oo' ( 79 ) 

M 2 

which now has the same form as at LL order. 



14 This formalism was first introduced at fixed order in 1351 
^99 



15 Notc that Agg has no small x poles at fixed order until NNLO. 
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Figure 8: Diagram contributing to the FF coefficient Cjf which is not included in the exact kinematics 
approximation. 



One must also consider the heavy flavour contributions to the longitudinal structure function. Equivalence 
of the VF and FF descriptions demands (in double Mellin space): 



C^ LL ( 7 ,N,M 2 /Q 2 )- 
C£ L H ( 1 ,N,M 2 /Q 2 )q H+ 



c FF,NLl s , v , ,-J, vl ;2 



\j,N,M 2 /Q 2 ) 



CV^h, N, M 2 /Q 2 ) + C v L F NLL h, N, M 2 /Q 2 ) 



Then use of the relations (|74I75[) yields: 



n VF,NLL _ n FF,NLL n LL aLL n VF,LL aLL 
°i,S _ U L,g ^L,H A Hg ~ ^L,g A gg ■ 



(80) 



(81) 



The term involving Agg must be neglected for the reasons given previously. Furthermore, the term in Ah 9 
presents a problem in the exact kinematics approximation as Q 2 /M 2 — > oo. This term diverges as (Q 2 /M 2 ) 7 
in double Mellin space, whereas the FF coefficient (coming from the exact kinematics massive impact factor) 
is finite as Q 2 — > oo. One can see this by noting that divergences in the longitudinal impact factor come 
from diagrams such as that shown in figure [8] Radiation of a gluon from the heavy quark leads to a collinear 
logarithm of type log(M 2 /Q 2 ). Such a diagram is a correction to the longitudinal impact factor that is not 
included by the exact kinematics approximation, and is needed to cancel the similar divergence coming from 
the term in An g - The upshot of this is that the right-hand side of equation (f8Tj) is divergent, whereas the 
left-hand side cannot be as one must have: 



lim 

Q 2 



c 



VF.NLL 
L,g 



r<NLL 
U L,g > 



(82) 



where the right-hand side denotes the massless longitudinal coefficient. Thus given the approximated impact 
factor does not have the correct asymptotic behaviour to match the term in Ajjg, this latter term must also 
be neglected along with the term in A gg . The heavy flavour contributions to Fl are in any case very small 
until higher Q 2 , where they approach the known massless terms. 



3.1 Approximate NLL Coefficients 

The results of the previous section can in principle be applied to find the NLL resummed heavy flavour 
coefficient functions. One uses the impact factors with exact kinematics to estimate the true NLL impact 
factors. However, as shown in [5D], it does not seem possible to obtain a closed analytic form for the 
exact kinematics impact factors in the case of heavy quarks. Power series results in N and 7 can be 
found numerically for given values of M 2 /Q 2 and used to investigate the NLL resummation. A problem 
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occurs in using the approach described in this paper, where equations (f3T)ll53"|) are used without any further 
truncations. The power series expressions for the coefficient functions are ill behaved up to rather high t, 
and parameterised numerical results cannot be obtained without a full definition of the NLL impact factor 
(i.e. not just the expansion in 7). So in this section we investigate the coefficient functions obtained using 
a gluon and structure function truncated separately to NLL order as described in [TU], but where the full 
gluon is used in the denominator when defining the coefficient function via: 

= (83) 

The power suppressed correction for these coefficient functions is much smaller at t values of phenomenolog- 
ical interest, and the conclusions drawn about how to approximate the resummed heavy flavour coefficients 
should be independent of the choice of NNLL terms. 

In [5D] it was noted that a smaller suppression of the heavy coefficient functions is induced by the NLL 
impact factor corrections than in the massless case. It is useful, in both the massless and massive sectors, 
to compare this suppression with that induced by the rest of the NLL framework. Overall, there are four 
factors which act to reduce the resummed splitting and coefficient functions at small x: 

• Running Coupling Corrections. 

• The implementation of the NLL BFKL kernel. 

• The use of the NLL gluon density Q\ instead of the LL result. 

• NLL contributions to the impact factors from exact kinematics. 

The effect of the running coupling in suppressing small x divergence has already been noted. It is also 
important to realise that the second and third mechanisms listed above account for most of the further 
suppression at NLL order, and the impact factors give a smaller correction once the kernel and gluon have 
been accounted for. 



One can see this in the massless case by looking at figure [3J which shows the LL result for P oa with running 
coupling corrections together with NLL results with and without the impact factor correction 16 !. At the lower 
t value, the effect of the NLL term in the impact factor is small compared with the combined suppression 
resulting from the gluon and kernel. At the higher t value, the impact factor induced suppression is relatively 
larger, but still reasonably small. Note that such a t value is relevant to the massive case, where one is only 
interested in the heavy flavour coefficients for t < 7. The VF coefficients decrease rapidly as t increases, as 
shown in [5] . 

One can examine the massive case in more detail using the impact factor results of 20J . It is not possible, 
however, to evaluate a power-suppressed correction once the exact kinematics are included as complete closed 
forms for the approximated impact factors are not known. Thus in order to investigate suppression effects 
one must choose a reasonably high value of t and hope that the correction is small. We consider the fixed 
flavour coefficient C^g, as the variable flavour coefficient will be tiny at values of t high enough for the 
power-suppressed correction to be small. The fixed flavour coefficient is shown in figure [TO] for t — 7. The 
LL power series result with running coupling corrections is again significantly bigger than the NLL result 
(although not as much as in the massless case), and the correction from the impact factor is much smaller 
than that due to the change in kernel and gluon density. The exact kinematics give a relatively larger effect 
at high x, but the fixed order NLO coefficient will dominate here. 

A similar result for Cf^ is shown in figure [IIJ Again the change induced by the NLL impact factor is of 

16 Note that the splitting function in these plots arises from a truncated gluon and structure function for comparison with the 
massive results. 



23 




io" 6 io" 5 io" 4 io" 3 io" 2 io _1 i io'So^io^io^io^io" 1 1 

X x 



Figure 9: The LL splitting function P qg for n/ = 4 with running coupling corrections (solid). Also shown is 
the NLL result with no impact factor correction (dashed) and the NLL resummed result which includes this 
correction (dotted). Note that the NLO splitting function has been added to give the correct behaviour at 
high x. 
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Figure 10: Power series results for the fixed flavour coefficient C^g for rif = 3 and t = 7. The LL result with 
running coupling corrections is shown (solid), alongside the NLL results with (dotted) and without (dashed) 
impact factor corrections. 
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Figure 11: Power series results for the fixed flavour coefficient C£ g for n/ = 3 and t — 7. The LL result with 
running coupling corrections is shown (solid), alongside the NLL results with (dotted) and without (dashed) 
impact factor corrections. 



a size that is small when compared to the suppression induced by the gluon and kernel. Thus, in approxi- 
mating the heavy flavour coefficient functions, it seems reasonable to ignore the contribution from the exact 
kinematics impact factors and only include the NLL kernel in the structure function and gluon density. Then 
one can parameterise the FF and VF coefficients over the whole desired range of t. Figure [9] suggests that 
the effect of the exact kinematics is felt more at higher values of t. This is encouraging for the approximation 
proposed here, as the variable flavour coefficients for decay to zero quite rapidly (~ M JQ ) as t increases. 

The results in this section have been derived using the NLL truncated structure function and gluon density. 
One can, however, check in the massless case that the suppression from the NLL impact factor contribution 
is also relatively small when the full structure function and gluon density (as described in section 2) are used. 
Figure [T2l shows the resummed P qg at t = 7 and rif = 4, with and without the impact factor contribution. 
One again sees that the dominant suppression in going to NLL order comes from the kernel and gluon rather 
than the impact factor. Given also that this procedure for generating the splitting function differs from that 
used throughout this section only by NNLL terms, the conclusion reached about the approximation for the 
massive coefficients seems to be robust. Indeed, impact factor corrections have a lesser effect in the heavy 
flavour sector, as noted previously. 

With this approximation for the heavy flavour quantities, the FF coefficients can be accurately parameterised 
in TV-space by a function of form (|48|) with positive powers of t. A similar function is sufficient for C% J , but 
with an overall factor of exp(— t) to reflect the decay of the gluon coefficient as t — » oo. There is a further 
subtlety involved in the longitudinal coefficient C\ F gl as it must tend to the massless resummed coefficient 
as Q 2 — ► oo. This will not happen if the NLL impact factor contribution in C\ F is ignored. To ensure 
that the correct asymptotic limit is reached, it is sufficient to model the NLL impact factor contribution 
by (1 — M 2 / Q 2 ) 15 where h^/ is the massless O(N) contribution to the impact factor, and the prefactor 
ensures a smooth matching between the fixed flavour coefficient at Q 2 = M 2 and the massless coefficient as 
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Figure 12: The LL splitting function P qg with running coupling for t = 7 and n/ = 4 (solid) compared with 
the NLL result with (dotted) and without (dashed) impact factor corrections, derived using the procedure 
outlined in section 2. 
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Figure 13: The approximate NLL resummed heavy matrix element An g for n/ = 4 and Q 2 = M 2 = 
2.25GeV 2 . 



Q 2 — > oo. Resummed heavy flavour coefficients for the quark singlet distribution for both F^ 1 and F^ are 
again approximated by using the colour charge relation (|60p . 

3.2 NLL Heavy Matrix Element A Hg 

The final ingredient for an approximate NLL analysis is the resummed matrix element governing the initial 
value of heavy quark distributions, evaluated at the matching scale Q 2 = M 2 . This is given in double Mellin 
space by equation (|78[) and one could use the asymptotic form of the impact factor, known in power series 
form from [20] . to obtain an analytic x space result according to equation (|53p . However, at LL order the 
analytic result was found to suffer from a large power suppressed correction [9]. In this paper we have also 
seen that NLL resummed quantities tend to have larger power suppressed corrections than their LL counter- 
parts, persisting to larger values of t. Hence we choose to define the heavy matrix element via the x-space 
equivalent of equation (|79|) corresponding to the approximate choice for the variable flavour coefficient. This 
is equivalent to evaluating A^ L by using the LL impact factor only with NLL kernel and gluon density 
corrections as has been done for the coefficients, and ensures continuity of the structure function across the 
matching scale in passing between the FF and VF descriptions. 

A plot of the heavy matrix element at Q 2 = M 2 for M = 1.5 GeV is shown in figure [T31 and one sees that it 
is negative over most of the range of x. This is in contrast to the LL result with running coupling corrections 
[S], which is positive as x -> 0. However, this result may or may not hold true once the complete NLL impact 
factor corrections (unknown) are included. Even so, the matrix element gives a very small correction to the 
structure function. 



27 



This now completes the set of quantities needed to perform an approximate NLL resummed global fit 
analogous to the LL fit of [9]. Further details of this fit and the results obtained are presented in the next 
section. 



4 A Global Fit at NLL Order 

The results of the previous section allow for a NLL+NLO analysis of DIS and related data. In order to 
compare with the LL+LO fit results of [9 , we use the same data sets as in that paper. We include F 2 
structure function data from the HI [TJGSIGIT] and ZEUS 2, 38, 39J collaborations at HERA; proton data 
from BCDMS [40], NMC [41], SLAC [H[43] and E665 [44]; deuterium data from BCDMS [45], NMC, SLAC 
and E665; CCFR data on F^ )N \x, Q 2 ) 06l[47]; data on the deuterium-proton ratio F 2 D /F£ from NMC [48]; 
charged current data from HI [T] and ZEUS [35]; data on the charm structure function F 2}C from HI [50] and 
ZEUS [EE]. The non-DIS data sets used are Drell-Yan (DY) data from the E866/NuSea collaboration [52] ; 
DY asymmetry data from NA51 [53]; data on the DY ratio cr^/c^y from E866 [54] ; W-asymmetry data 
from CDF [55]. We also check our predictions against the Tevatron jet data [Ml [57], although these are not 
explicitly included in the fit. In order to compare with the fixed order QCD expansion, we also produce a 
fit at NLO in the DIS scheme. Before presenting results, we first discuss some implementation issues arising 
in the fitting procedure. 



4.1 Implementation 

Fixed order cross-sections for each of the data sets mentioned above must be transformed into the DIS 
scheme. For most quantities, the transformation of splitting and coefficient functions can be carried out 
analytically. For others, some parameterisation is necessary. For example, the NLO fixed flavour (FF) 
neutral current coefficient functions for the heavy contribution to the structure functions transform as: 

CZ mm = C^^« + 2n<f «) ® cgM ; (84) 

r ,FF,DIS{2) _ r ,FF,MS(2) , r rFF,MS(l) « ^MS(l) fnc\ 
^ 'a.q a.q ' ^a,g w a.q J \°'- ) / 

where a € {2,L} and the rif dependence has been explicitly displayed in the light gluon coefficient. The 
convolution on the right-hand side cannot be calculated analytically. Furthermore, the NLO MS scheme 
heavy flavour coefficients are also not known analytically [58, 59 , although recently a new parameterisation 
has been produced 60J. We parameterise the transformation terms in equations (|84I85[) by evaluating the 
convolution in each case numerically for a range of values of x and M 2 /Q 2 and fitting these to a function of 
form: 



^ I Q 2 



n . 111 



m 2 \ r r . , m ( x 



a n +b n (x- X m ax)] log™ + C n (x - X max ) m > , (86) 



where x max = (1 + AM 2 / Q 2 )^ 1 . This functional form then satisfies the kinematic constraint coming from 
the fact that the partonic centre of mass energy must exceed the threshold for heavy quark pair production. 
Fortran code for the DIS scheme NLO heavy flavour coefficients is available on request. We define the vari- 
able flavour (VF) heavy flavour coefficients according to the DIS(x) scheme [9j. 

The charged current heavy flavour coefficients are more complicated, as there are 0(a° s ) contributions to 
the structure functions resulting from the production of a single charm quark in the final state: 

F 2 C = 2[cos 2 6> c £s(£) + sin 2 c £<i(O] = C^J ® s + £7$ g d . ( 87 ) 

xFI = 2[cos 2 6 c xs{£) + sin 2 9 c xd(£)] = x[C { 2 °] ® s + dfj ® d] 

where C is the Cabibbo mixing angle and £ = x/x^ x . The upper limit: 

Q 



<ax = [ 1 + 7>f ) ( 89 ) 
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is the charged current analogue of the variable x max introduced above for neutral current scattering. We here 
present results for the transformation of the NLO FF coefficients for charged current scattering, as they have 
not been given in the literature before. First one must carefully define what is meant by NLO in this case. 
Denoting the weak eigenstate [cos 2 9 c £s + sin 2 9 c ^d] by s, the production of s (the weak eigenstate conjugate 
to c) vanishes at zeroth order for W + exchange. What is conventionally regarded as NLO in the fixed order 
expansion includes NLO production of this conjugate eigenstate, and thus includes the 0(a 2 s ) FF coefficient 
functions. Then results up to this order for the transformation of the FF coefficients are (i £ {2, 3}): 

C FJ(X)m S = c iJ(l)MS _ C (1)MS ^ (c (0) + C (0) ); (gl) 

C<£f (2)DIS = C<£f (2)M5 + 2njCfJ [im ® Cgf 5 - ® (G< j + C!$) 

- 2nAT ® ® iCg + C$) - ® (C(f is + cg DIS ) 02) 

C FF(2)DIS = C FF {2) MS + c mi)m C (1)MS _ J_ C g)MS _ ^MS g ^(1)5*3 g (c (0) + ^ (g3) 

n/ 

Analogous results apply in the case of M^" exchange, where production of s vanishes at LO. Some care is 
needed to obtain the above results, as one must use the MS — * DIS transformation of a single quark species, 
which is: 

q? 1S = C™ s ® <zF + JLtfgJ fi) + ® s™, (94) 

with a similar equation for antiquarks. Some simplification of the above results is possible. In [61 it is 
noted that the NLO correction to single charm production, represented by C^J^, is negligible in practice, 
and can be approximated by the massless quark coefficient. In the DIS scheme, this is zero and so can be 
neglected. Note this also means that one can neglect the final term in equation 192(1 . 

A slight problem arises in that the NLO quantities (jFF(2)US are not known for charged current scattering. 
In [60] they are approximated by evaluating the neutral current results with modified threshold behaviour 
i.e. dependence on x max replaced with dependence on x^ x . We follow this approach here, although the 
transformation terms to the DIS scheme can be parameterised using exact results given that the LO quan- 
tities are all known. 



The variable flavour coefficients are only needed up to O(as), and adopting the DIS(x) scheme choice as 
outlined in [9] gives: 

„VF(1)DIS _ n (i)MS ^(l)MS ^(0) . / QC .s 

C VJ(1)DIS = C VJWm _ ^(0) ^ C (1)MS _ C (1)MS ^ (c (0) + C (0) } _ (g6) 

Approximation of the NLO single charm coefficient by its massless counterpart means that one can neglect 
the contribution from equation (|95p in the DIS scheme. 

One also needs to specify a scheme for the NLO coupling when heavy flavours are involved. We use the 
prescription of Marciano [62j . The NLL approximate resummed splitting and coefficient functions have been 
determined by solving the BFKL equation at NLL order but with the LO running coupling. It is not then 
correct to use the NLO running coupling in these expression^, and in the fitting program one must define 
the LO coupling by matching with the NLO coupling at some scale. We use Q 2 — 16GeV 2 , as this is 

17 The overall normalisation in the impact factors does contain a power of the NLO coupling, as the scale of this coupling is 
Q 2 and not k 2 when solving the BFKL equation (see |10|). 
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Figure 14: The function (1 — x) 30 used to filter out high x information in the alternative NLL resummed fit. 



guaranteed to correspond to four active quark flavours. One then has: 



I _ ^NLO {t NL 0) (97) 



[3ot LO ' 

which defines the LO t parameter at this scale, where as usual t — \og(Q 2 / A 2 ) 



For the DY data it is conventional to implement the NLO corrections to the cross-section via a multiplica- 
tive if-factor. Rather than produce a if- factor in the DIS scheme, it is sufficient to explicitly transform the 
parton densities to the DIS scheme in the existing MS scheme code. We adopt such an approach here. 

As noted in section \TM some of the NLL improved splitting and coefficient functions contain terms at O(N ) 
and O(N) in Mellin space, which we model by a suitable function of x in practice. It is perhaps unnatural 
that the small x resummation should lead to such noticeable effects at large x. Thus, for comparison with 
the standard NLL resummed fit, an alternative fit is also presented in this chapter in which the all resummed 
quantities have been multiplied by a factoiF^I (1 — a;) 30 . A plot of this function is shown in figure HU and 
one can sec that its effect is to filter out high x information in the resummed quantities for x > 10~ 2 . 



4.2 Results 

The x 2 values for each data set obtained in the NLL and alternative NLL fits are shown in table [1] alongside 
results obtained using a NLO DIS scheme fit with no resummations. We also include for comparison the LL 

18 One may choose any power of (1 — x) for this damping function that is sufficient to filter out the high x behaviour. Note 
that a function of form (1 — x) leads to a cutoff of the resummation for x > N -1 . We have checked that N = 20 gives very 
similar results. 
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Data Set 


No. data pts 


Xnll 


X.NLL(2) 


Xnlo,dis 


xi L 


HI ep 


417 


352 


413 


415 


341 


ZEUS ep 


356 


273 


245 


239 


273 


F% 


27 


25 


43 


40 


25 


BCDMS ft p 


167 


171 


158 


180 


186 


BCDMS (jt D 


155 


198 


202 


219 


223 


NMC \i p 


126 


126 


131 


139 


120 


NMC \i D 


126 


102 


103 


103 


97 


SLAC ep 


53 


65 


78 


76 


92 


SLAC eD 


54 


53 


81 


65 


82 


E665 /i p 


53 


61 


60 


61 


57 


E665 fi D 


53 


51 


58 


66 


54 


CCFR F% N 


74 


74 


75 


62 


97 


CCFR Fg N 


105 


115 


104 


114 


138 


HI CC 


28 


29 


31 


32 


32 


ZEUS CC 


30 


44 


47 


49 


44 


NMC n/p 


156 


147 


165 


154 


157 


E866/ NuSea DY 


174 


240 


225 


234 


296 


NA51 DY asym. 


1 


5 


5 


13 


16 


E866 o&/o% Y 


15 


7 


7 


14 


7 


CDF W asym. 


11 


16 


16 


17 


20 


Total 


2181 


2154 


2247 


2289 


2357 



Table 1: The quality of fit from the NLL and modified ((1 - a;) 30 weighted, labelled NLL(2)) fits for each 
dataset, as well as results from a NLO DIS scheme fit and the LL fit of [9]. 



fit values described in [9] . The numerical results quoted in that paper are erroneous due to a mistake in the 
fitting code. The LL values in table [T] are obtained with the corrected program. A feature of the previous 
LL fit, not explicitly noted in [9], is the need for a significant violation of momentum conservation due to the 
fact that the LL resummcd splitting functions do not satisfy the momentum sum rule. The input partons in 
the LL fit of table [1] carry a total of 78% momentum, consistent with the fact that the evolution is enhanced 
over a large range of x. All of the next-to-leading fits in table [IJ on the other hand, are obtained with full 
momentum conservation imposed. The input partons carry 100% momentum, and we saw earlier that the 
NLL resummed splitting functions satisfy the momentum sum rule to a very good approximation. 

All three next-to-leading fits are good all-round fits. Of the two resummed fits, the standard NLL fit clearly 
outperforms the NLO DIS-scheme fit, whereas the fit using the (1 — a;) 30 weighting gives less of an improve- 
ment. The fit to the HI data in particular suffers without the high x part of the resummation, indicating 
that this is indeed needed for a truly good fit. Comparing the NLL resummed fit with the NLO results, it 
seems that the introduction of rcsummations seem to decrease the tension between different data sets - even 
though one might expect some tension from not including resummations in the coefficient functions for all 
observables. The NLL fit is much better than the LL fit. One sees a more satisfactory description of the DY 
data (due to the NLO corrections), and also the fit to the deuterium data improves. This is a hint that the 
sea and valence distributions have the correct relative behaviour, and thus that the suppression of small x 
resummations anticipated in [9] has indeed occurred. 

The F2 values from the fit at small x are shown with the data in figures [15] and HH As already noted 
at LL order [9], one sees that resummations are needed to give the desired slope in F2 as Q 2 increases, as 
is particularly evident in figure 1161 Curiously, the modified NLL fit with the high x part of the resum- 
mation filtered out (labelled NLL (2) in the plots and tables) gives a very similar prediction for F2 to the 
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Figure 16: Theoretical predictions for the structure function F 2 alongside the data, for 5 x 10 4 < x < 
8 x 1CT 3 . 
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Fit 




NLL 
NLL(2) 
NLO 


0.1184 
0.1130 
0.1181 



Table 2: The values of the strong coupling resulting from each of the next to leading fits. The world average 
(at NLO) is 0.1187(20) [53]. 



NLO fit. One expects the theory prediction at 10~ 4 < x < 10~ 3 to be sensitive to the coefficient and 
splitting functions at 10 -2 < x < 10" 1 through the convolutions, and so it is perhaps no surprise that 
filtering out the resummations in this latter regime greatly diminishes their effect at lower x. The modified 
fit only converges with the full NLL fit at very low x ~ 10 -5 . This suggests why the modified fit is less 
successful than the full NLL fit - it differs slightly from the NLO approach in the region where the fixed order 
fit works well, and does not contain enough of the small x resummation to improve matters very well at low x. 

The charm data is shown in figure [T7J Once again, the modified resummed and NLO fits give very similar 
predictions until one reaches very low x. The full NLL lies above the fixed order prediction only at the 
lowest x values, and otherwise wants to tend towards the lower limit allowed by the data. Nevertheless, the 
NLL fit to the charm data is excellent as can be seen from tabic [T] 

The strong couplings as(M§) required by each fit are given in table O The NLO and full NLL results are 
compatible with the world average, whereas the NLL(2) coupling is rather lowr^l. The latter result is not 
very precisely determined, however. For example, a value as(Af|) = 0.116 gives a very similar fit (with 
minor reshuffling of the gluon parameters) in the NLL(2) program. That the modified fit gives the lowest 
value is expected given that the modified resummed splitting functions have an enhancement at low x, but 
have less of a dip at moderate x as this has been filtered out (see figure [3j . Thus a lower coupling is needed 
to fit the data at small x. 

The improved fit to scattering data suggests that in the NLL fit the resummations are confined to lower 
values of x. To investigate this in more detail, one may examine the gluons obtained in each of the fits. 
These are shown in figure QS for Q 2 = lGeV 2 and Q 2 = 100 GeV 2 . At the input scale Q 2 = lGeV 2 the fixed 
order DIS-scheme gluon, although positive at very small x, is negative over a large range of x. On the other 
hand the resummed gluons are positive in this regime, with the full NLL gluon rising at small x. Although 
positivity of the gluon is not necessary, such a feature is somewhat more physically appealing (for example, 
structure functions will then be positive at low Q 2 as long as the coefficient functions are). At moderately 
high x, the modified NLL gluon is very similar to the NLO result, as expected from filtering out the small 
x resummations here. The gluon from the modified fit also lies below the full NLL gluon at small x, which 
is to be expected given the increased dip at moderate x in the full NLL splitting functions. 

Looking at the right-hand panel of figure [HI one sees that at higher Q 2 the three gluons are extremely 
similar at moderate and high x, with a noticeable difference occurring only for x < 10~ 2 . The full resummed 
gluon lies above the fixed order result as x — > due to the increased small x evolution, and the gluon from 
the modified fit is closer to the fixed order gluon than the full resummed gluon. This is consistent with the 
above noted fact that the predictions for F2 in the fixed order and modified resummed fits are very similar 
until very small x. 

Despite the good fits to the DIS and related data from all of the approaches considered here (with the full 
NLL fit giving a global y 2 of less than one per point), there is a worrying feature in the input gluon at high 
x. Both the fixed order and resummed gluons show a negative dip at very high x. This behaviour of the 

19 Strictly speaking, the couplings in the resummed and fixed order fits are not the same quantity. 
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Figure 17: Resummed predictions for the charm structure function F% alongside HERA data, for 1.3 x 
i<2x 1(T 2 . 
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NLL resummation with negative high x gluon and DO jet data, % = 303/82 pts 



NLL resummation with negative high x gluon and CDF1B jet data 
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Figure 19: The comparison of theory to the jet data for the NLL parton distributions using coefficient 
functions at NLO in as- For the DO data [S5] (left) the blue points represent the size of the systematic 
errors while the uncorrelated errors are shown explicitly. For the CDF data [57] (right) the movement of 
theory relative to data using the systematic errors is shown. 
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gluon is not consistent with the Tevatron jet data, which prefers at least a positive semi-definite gluon at 
high x. It was shown in [64] that a positive DIS-scheme gluon at high x automatically gives a good fit to 
the Tevatron jet data. We now see that a negative gluon certainly does not. The comparison to the jet 
data achieved by transformation of the gluons from the above fits to the MS scheme gives Xjets — 300 or 
worse for 113 data points. The quality of the comparison is shown for the NLL partons in figure fT9l The 
negative high- a; gluon in the DIS scheme is in fact still slightly negative in the MS scheme, and as such 
at high values of Et and/or high rapidity the theoretical prediction is very much below the data. This is 
similar, but slightly better, for the NLL(2) and DIS-scheme NLO partons in this section. In the case of the 
DO data the resulting x 2 for all the parton sets is very poor indeed, i.e. \ 2 ~ 300 for 82 points. For the 
CDF data it is not so much worse than the standard good fits \ 2 ~ 70 rather than 50 for 31 points, but 
this is achieved only by assuming that the systematic errors are conspiring to make the true data very much 
smaller than the central value presented. Clearly this type of high-x gluon distribution in the DIS scheme is 
not acceptable. The solution to this problem is considered in more detail in the next section. 

4.3 Global Fit with Positive Large x Gluon 

In the NLO fixed order fit one can constrain the gluon input parameters to ensure a positive high x gluon. 
Subsequent refits will converge to the local minimum satisfying this condition. In the full NLL resummed 
fit, fixing the gluon parameters is not sufficient to give the desired effect in the gluon, due to the term 
proportional to 5(1 — x) in P gg which is negative in practice. By the colour charge relation of equation 
(|60|) . this gives a negative delta function contribution to P gq leading to a negative evolution in the high x 
gluon whose derivative at a given Q 2 becomes proportional to the quark singlet. This is a relatively large 
effect given that the gluon dies away at large x much more quickly than the quarks. However, there are 
strictly speaking no 5(1 — x) terms in the full fixed order P gq - the resummed delta function is an artifact 
of expanding about N — 0. One can replace the delta function term with a function that has the same first 
moment and little effect at small x. A suitable choice is jV -1 , as the coefficient of this term is sufficiently 
small to have negligible effect at low x. 

We consider the full NLL resummed fit together with a fit with the resummations weighted by a (1 — x) 30 
factor as before. The resulting \ 2 values are shown in table [3] where the fits are labelled with a plus sign 
to distinguish from the previously described next to leading results. Both the fixed order and full NLL fits 
are very good overall. The resummed fit performs better, with the main improvement occurring in the HI 
(including charged current scattering) and charm data. It is not possible to fit the HI data as well as in the 
previous NLL fit once a positive gluon is required - the HI data would like a lower gluon at high x, at the 
expense of a good fit to the Tevatron jet data 20 !. The improvement in the NLL fit over the fixed order result 
is considerable, with an overall \ 2 difference of ~ 100 between the two fits. The modified resummed fit again 
shows much less improvement. The structure function F% from the fit is compared with the small x data in 
figures I20H221 Again the slope of the resummed prediction is higher at small x as is required by the data, 
but now the effect is less dramatic. The resummed prediction for the charmed data is again lower for most of 
the x range, but gives an excellent fit. The modified resummed results are too close to the NLO results until 
very small x, and thus cannot provide as good a description. One may also compare the gluons obtained 
from the fits described here (shown in figure [23)) with those of figure [TBI The gluons at Q 2 — 100GeV 2 are 
much the same as before (except for high x), and mutually similar for x > 10~ 2 . Pleasingly, the gluons in the 
left-hand panel of figure [23] are also very alike at high x suggesting that the small x resummations are indeed 
not affecting the high x behaviour of the gluon (hence, the structure functions). The NLO gluon is negative 
for x < 10~ 2 , turning over at very low x but not becoming positive over the range of interest. On the other 
hand, the full NLL resummed gluon is still positive at small x and growing as x — > 0, albeit very gently. 
Thus the resummed gluon is now positive definite over the whole range of x. The gluon from the modified 
fit is, as expected, somewhere between the full resummed and fixed order gluons and is now negative (but 
tending to zero) as x — > 0. This is a further indication of the inadequacy of the modified approach - it now 

20 This has been noticed previously by global fits at NLO e.g. [6511661 . 
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Data bet 


No. data pts 


■2 

X.NLL+ 


2 

%NLL(2) + 


2 

X.NLO.DIS+ 


TT1 ~ 

HI ep 


417 


oca 

360 


A O O 

438 


A OO 

423 


ZEUS ep 


356 


260 


zol 


O A O 

242 




07 


ZD 


6Z 


4U 


BCDMS m P 


167 


166 


188 


186 


BCDMS /U D 


155 


228 


217 


231 


NMC fi p 


126 


126 


140 


130 


NMC fi D 


126 


95 


106 


100 


SLAC ep 


53 


74 


72 


73 


SLAC eD 


54 


58 


75 


64 


E665 n p 


53 


62 


59 


59 


E665 fj, D 


53 


49 


58 


61 


CCFR ^ 
CCFR Fg N 


74 


58 


62 


67 


105 


160 


124 


153 


HI CC 


28 


30 


31 


31 


ZEUS CC 


30 


45 


44 


45 


NMC n/p 


156 


163 


153 


167 


E866/ NuSea DY 


174 


255 


241 


246 


NA51 DY asym. 


1 


10 


5 


7 


E866 o v ^ Y la^ Y 


15 


9 


7 


9 


CDF W asym. 


11 


15 


17 


18 


Total 


2181 


2249 


2318 


2352 



Table 3: The fit quality from the full NLL resummcd, modified rcsummed and fixed order DIS-schcme fits, 
where the gluon is constrained to be positive at high x as is required by the Tevatron jet data. 
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Q 2 (GeV 2 ) 



Figure 20: Theoretical predictions for the structure function F 2 alongside the data, for 5 x 10 5 < x < 
4 x 10~ 4 , with the gluon constrained to be positive at high x. 
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Figure 21: Theoretical predictions for the structure function F 2 alongside the data, for 5 x 10 4 < x < 
8 x 10~ 3 , with the gluon constrained to be positive at high x. 
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Figure 22: Resummed predictions for the charm structure function F% alongside HERA data, for 1.3 x 
i<2x 10~ 2 , where the gluon is constrained to be positive at high x 
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NLL resummalion and DO jet data, a s (M z )=0.119 , x 2 = 68/82 pts 



NLL resummation and CDF1B jet data 
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Figure 24: The comparison of theory to the jet data for the NLL+ parton distributions using coefficient 
functions at NLO in as. 



agrees closely with the NLO framework at high x, but has none of the benefits of the small x resummation 
at moderately low x. That the gluons are now qualitatively different in the full and modified resummed fits 
indicates that the former approach is more appropriate. Note that the gluons in the right-hand panel deviate 
from each other at x ~ 0.005. This is consistent with a previous analysis which added phenomenological 
log ra (l/a;) terms at C(a s ) to the splitting functions, and found an input gluon qualitatively similar to 
that obtained by imposing a cut x > 0.005 on the data, which thus gives some indication of the value of x 
at which resummations become important. 

The fit quality for the jet data is now very good, as expected from the gluon modification. After transforming 
the gluons to the MS scheme, all three fits give x 2 — 120 for 113 points (about as good as is possible in a 
global fit) with the full NLL fit giving xjets = H9. The comparison of the predictions using these partons to 
the data is shown in[Ml and gives x 2 = 68/82 for the DO data and \ 2 — 51/31 for the CDF data. The situ- 
ation is very similar for the other two parton sets. We note that the quality of the fit to the rest of the data 
has deteriorated by 50 — 100 units of \ 2 by the constraint that the high- a; DIS scheme gluon be positive, but 
the improvement in the fit to the Tevatron jets data is over 200 units better, so it is clear that the imposition 
of the positive high- a; DIS scheme gluon is demanded by a global fit. The fact that the jet data has not been 
explicitly implemented in the fit is further evidence that a positive semi-definite gluon in the DIS scheme 
automatically leads to a good fit to the Tevatron data in the MS scheme [64]. Indeed, the gluon in the DIS 
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Fit 


as (MI) 


NLL+ 
NLL(2)+ 
NLO+ 


0.1188 
0.1196 
0.1195 



Table 4: The values of the strong coupling resulting from each of the next to leading fits with the high x 
gluon constrained to be positive. The world average (at NLO) is 0.1187(20) 63J. 



scheme is still very small at high x, i.e. ~ (1 — x) 71 , where 77 = 7 — 10 in the various sets. It is then much larger 
in the MS scheme. Actually including the Tevatron jet data in the fit can lead to an further improvement 
of 5 — 10 units for this data, with a slightly larger high- a; gluon. However, this is at the expense of a de- 
terioration of 5 — 10 units for the rest of the data, so the quality of the global fit is much the same in each case. 

The values of as obtained in the positive gluon fits are shown in table |H The NLO coupling has increased, 
which makes sense given the smaller gluon at moderate and small x. The full resummed fit result is in 
excellent agreement with the NLO world average, despite the tendency of resummations to want to decrease 
the coupling. The value of ag(M|) from the modified fit is very similar to that obtained in the NLO fit, as 
perhaps is expected from filtering out resummations in the moderate and high x regimes. 



5 NLL Resummed Prediction for Fl 

With gluons that are more consistent with the Tevatron jet data, we are now in a position to examine the 
longitudinal structure function. The predictions for Fl using the three global fits of the previous section are 
shown in figure 1^51 The NLO DIS-scheme result, unlike the corresponding MS-scheme result [8], is positive 
at low Q 2 as x — > 0. However, it still dips slightly negative at moderate x which is unphysical. The NLL 
resummed result is much more sensible at low x, and lies above the NLO result, as can be expected from the 
increased gluon and coefficient function. The modified resummed approach leads to a prediction for Fl that 
is very close to the full resummation value at high x as expected, but lower at small x due to the smaller 
gluon. It catches up with the full resummed result as Q 2 increases due to the raised coupling, and shares 
with the NLO prediction the unsatisfactory feature of being oscillatory (though not nearly as much) at small 
Q 2 due to the negativity of the gluon. This is no longer a problem in the full resummed fit, which predicts 
that Fl is rising smoothly as x — > 0. 

We also show in figure [25] the LL result (with running coupling corrections) of [9] , corrected from that 
paper due to a mistake in the evolution code. One sees that the LL prediction is turning over as x — ► 0, 
becoming negative at low enough x and Q 2 . This is due to the negativity of the gluon in the LL fit, cured 
by NLL resummation in the present paper. The LL result dramatically undershoots the fixed order results 
at moderately high x, due to the prevalence of LL small x resummations in that region and to quite large 
contributions from the NLO coefficient functions. 

One may also compare the resummed prediction with the MS scheme fixed order predictions for Fl- Al- 
though Fl is formally scheme independent only when considering all orders in the perturbation expansion, 
comparing the resummed DIS-scheme prediction with the results obtained using the first three orders of the 
MS scheme expansion gives some indication of the stability of the resummed result. This is shown in figure 
l26l The resummed result lies above the MS scheme NLO and NNLO results at intermediate Q 2 , again due 
to the increased gluon in the moderate x regime (NLO MS gluon is negative at moderate and small x - see 
[63]). At higher Q 2 it is similar in magnitude to the MS scheme NLO result and one can clearly see that 
the problem of severely undershooting the fixed order results at moderate x has been alleviated by going to 
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Figure 25: NLL resummed predictions for the longitudinal structure function compared with the NLO 
DlS-scheme fit. 
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Figure 26: NLL resummed prediction for the longitudinal structure function compared with the fixed order 
MS scheme results. 
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Figure 27: NLL resummed prediction for the longitudinal structure function compared with fixed order MS 
scheme results. 



NLL order in the resummatiorl^l. 

The behaviour as Q 2 increases at fixed x is shown in figure [271 It lies somewhere between the MS scheme 
NLO and NNLO results at very low Q 2 and x, and is somewhat flatter as Q 2 increases due to the dip in the 
gluon splitting function at intermediate x, to the reduced coupling in the resummed fit, and to the fall off 
of the resummed coefficient function. 



6 Discussion and Conclusions 

In this paper the approximate framework for NLL small x resummations introduced in [HO HO] has been 
used with the approach of [TUl [9] to produce NLL+NLO splitting and coefficient functions in both the mass- 
less and heavy sectors. Our results for the splitting function P gg are very similar to those of alternative 
approaches [211 E2] , even though we do not resum collinearly enhanced terms in the BFKL kernel [24] . We 
stress, however, that in processes where there are two perturbative scales at either end of the BFKL 4-point 
function such a kernel resummation would indeed be appropriate. 

We have successfully applied our results to scattering data, and a new DIS-scheme global fit at NLO in the 
fixed order QCD expansion has been provided for comparison. An excellent overall fit is obtained with the 
resummations, and a marked improvement over the fixed order results is observed. Particularly appealing is 
the positive definite gluon obtained in the NLL fit. Whilst not a necessary condition for the gluon density, a 
positive gluon avoids the problem of negative structure functions - exhibited by the NLO prediction for Fl , 

21 The fact that the resummed Fl undershoots slightly the fixed order results in figure [26] is due to the different factorisation 
schemes, as can be seen from figure [25] where this effect is lessened considerably. 
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which in both the MS and DIS schemes is negative for some range of x. Indeed, even if the overall fit quality 
were the same for the NLL and NLO fits one would probably consider the NLL fit as the better descrip- 
tion given the more physical form of the gluon. This is also evident in the prediction for Fl, in which the 
NLO fit shows an oscillatory behaviour at moderate x whereas the resummed prediction grows more sensibly. 

This more physical gluon also provides a more sensible starting point for an investigation of nonlinear cor- 
rections to the gluon evolution at small x [67j [68] . It seems most appropriate to approach this problem 
considering additional corrections on top of the, in comparison, well understood leading twist framework 
where in most regions of parameters space the partons are determined with great precision. This has al- 
ways proved to be problematic when using as a starting point the NLO or NNLO gluon distributions from 
global fits since these always become very small, or frequently negative, at Q 2 ~ 2 GeV 2 and x ~ 10 -4 . 
Consequently saturation due to large parton densities is largely counterintuitive. Sometimes LO in as 
partons are used, but these are qualitatively inaccurate at small x and Q 2 . The small x gluon produced 
here represents a more complete picture of the leading twist partons at very small x and should be a more 
physical quantity to work with when looking at additional corrections, particularly since many of the recent 
approaches are extensions of the leading twist physics within the BFKL framework [Ml [701 ITTl IT2] . The 
gluon in this paper would also be more appropriate for use in approaches which need unintegrated gluon 
distributions, e.g. [73], rather than obtaining this from fixed order perturbative gluons. The latter have not 
been obtained from the global fit in the same manner in which they are being used, whereas the gluon in 
this paper has been obtained by performing kx integrations over impact factors in a more compatible manner. 

The best fits obtained to the data sets produce a negative gluon at high x. Such a feature is not consistent 
with the Tevatron jet data which wants a positive semi-definite gluon in the DIS scheme at high x, leading 
to a large gluon in the MS scheme. One can constrain the gluon at high x to be positive semi-definite, and 
a very good overall fit is still obtained. The comparison to the jet data is then excellent, even though the 
Tevatron data is not explicitly included in the global fit. Again, the presence of small x resummations is 
noticeably preferred. 

The resummed splitting and coefficient functions contain important high x terms, which are really artifacts 
of an expansion about N = 0. Hence it is tempting to filter out these large x terms with a function such 
as that shown in figure [T4l However, in each of the fits this approach does not perform well, giving only a 
slight improvement over the fixed order description. The results are similar to, but slightly different from, 
the NLO results at moderate and high x. Predictions for structure functions do not approach the full NLL 
results until very low x. Thus there is not enough benefit from the small x resummation. Also, in neglecting 
the Dirac-type terms from the resummed splitting functions, important information is left out. The HI data 
in particular benefits from inclusion of the high x terms in the resummed splitting functions. Therefore, it 
is not really correct to think of this information as being at "high x" . It really is an artifact of a resummed 
expansion about N = 0, and all of the resummation thus obtained must be implemented. 

The NLL resummed fit is also better than the LL resummed fit of 0. The global fit results, gluon distribu- 
tion and Fl prediction all indicate that the shortcoming of the LL approach - that the resummations were 
too prevalent at moderate and high x - has been eradicated. Indeed, the LL fit is somewhat misleading in 
that a large momentum violation (—22%) in the input partons was necessary. All of the fits described in this 
paper have input partons which carry 100% momentum, and the evolution in the NLL evolution conserves 
momentum to a very good approximation^!. 

The fit to the charm data is greatly improved by implementing the (NLL resummed) DIS(x) scheme. The 
resummed prediction for F^ lies towards the lower end of the range allowed by the data except for at very 

22 We have explicity checked momentum conservation by integrating the partons at various values of Q 2 . We find that the 
momentum sum rule is violated by no more than 3.5%. The effect is cumulative and this value is only reached for the highest 
Q 2 . For much of the small x HERA range the Q 2 is such that the violation is less than 1%. 
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low x, but leads to an excellent fit. In fact the features are similar to those seen when including the NNLO 
corrections [60] . 



Further improvements to the NLL resummation are possible. The massless splitting and coefficient functions 
can in principle be updated once the full NLL impact factors are known. However, we have seen that this 
is not the most important correction. Furthermore, the massive coefficient functions may benefit from im- 
proved modelling of the NLL impact factor effects (although these have a much smaller effect on the overall 
fit). Nevertheless, the conclusion that small x resummations are beneficial in a global fit to scattering data is 
a firm one, made even more so by the (slightly) approximate nature of the approach used here. The partons 
are available from the authors in grid form in the direct DIS scheme and also in the MS scheme obtained via 
a fixed order transformation. However, these would only be physically meaningful at large Q 2 or for a range 
of (high) x where resummation effects in the MS — > DIS transformation are expected to be less important. 

To conclude, we believe the form of the gluon together with the x 2 improvement in the resummed global fit 
provide extremely compelling evidence that small x resummations are a necessary and beneficial prerequisite 
for describing the HERA data, and for predicting high energy partons for the LHC. 
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